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I.   INTRODUCTION 

A.      BACKGROUND 

The  thermodynamic  growth  of  sea  ice  is  essentially  a  one-dimensional 
heat  conduction  problem  with  the  boundary  conditions  expressed  in  terms  of 
incident  heat  fluxes  of  shortwave  and  longwave  radiation  and  sensible  heat  at 
the  surface  and  the  upward,  sensible  heat  flux  from  the  ocean  (Untersteiner, 
1961).    The  heat  budget  is  approximately  one-dimensional  because  sea  ice  is 
very  thin  in  comparison  to  the  large  area  it  covers,  and  small  changes  in  the 
atmospheric  or  oceanic  heat  fluxes  can  cause  major  changes  in  sea  ice 
thickness  and  extent  (Maykut,  1985).  The  surface  temperature  of  the  ice  will 
vary,  according  to  Kirchoffs  Law  (Stewart,  1985),  so  the  outgoing  longwave 
radiation  nearly  balances  the  incoming  fluxes  at  the  surface.   When  the  surface 
temperature  of  the  ice  is  less  than  the  freezing  temperature  of  the  ice  a 
temperature  gradient  is  set  up  across  the  ice  and  heat  heat  will  be  conducted 
upwards  through  the  ice.    If  this  conductive  flux  exceeds  the  oceanic  heat  flux 
due  to  entrainment,  latent  heat  will  be  released  at  the  bottom  of  the  ice  in  the 
form  of  ice  growth  (Thomdike,  1992). 


Early  thermodynamic  sea-ice  models  (Untersteiner,  1964;  Maykut  and 
Untersteiner,  1971;  Washington  et  al.,  1976)  used  observed  or  calculated 
atmospheric  forcing  with  prescribed  upward  ocean  heat  flux  into  the  ice.  The 
next  development  was  coupled  ice-ocean  models  to  include  ocean  dynamics. 
Parkinson  and  Washington  (1979)  used  a  mixed  layer  of  fixed  depth  and  a  fixed 
value  for  the  upward  oceanic  heat  flux.   Hibler  and  Bryan  (1987)  coupled  a 
dynamic  thermodynamic  sea  ice  model  to  a  multi-level  baroclinic  ocean  model, 
and  Semtner  (1987)  coupled  a  dynamic  thermodynamic  sea  ice  model  to  a 
multi-level  primitive  equation  ocean  circulation  model.   Both  of  these  later 
models  calculated  the  upward  oceanic  heat  flux  by  comparing  the  temperature 
difference  between  the  ice  and  the  upper  level  of  the  ocean  model.  Although 
full  mixed  layer  dynamics  are  not  used  in  these  models  to  calculate  the  upward 
oceanic  heat  flux,  they  do  show  the  sensitivity  of  ice  thickness  and  extent  to  the 
ocean  heat  flux.     It  was  not  until  Lemke  (1987)  coupled  a  thermodynamic  sea 
ice  model  with  a  one  dimensional,  ocean  mixed  layer  model  that  the  vertical 
oceanic  heat  flux  was  determined  using  mixed  layer  physics.   However,  his 
model  empirically  paramaterizes  the  mechanical  mixing  due  to  ice  motion. 
Mellor  and  Kantha  (1989)  also  developed  a  coupled  ice-ocean  mixed  layer 
model  utilizing  second  order  closure  to  determine  ocean  eddy  conductivity. 

Claes  (1990)  showed  the  impact  of  mixed  layer  dynamics  on  the  onset  of 
freezing  and  determined  that,  under  certain  conditions,  the  upward  oceanic  heat 
flux  could  prevent  ice  formation  regardless  of  the  atmospheric  forcing. 


B.      OBJECTIVE 

The  fluxes  of  momentum,  mass,  and  heat  into  the  ocean  are  controlled  by 
the  oceanic  surface  mixed  layer  so  proper  understanding  of  mixed  layer 
dynamics  is  critical  to  understanding  the  physical  oceanography  of  the  arctic. 

For  this  research  a  simplified,  sea  ice  model  (Thomdike,  1992)  is  coupled 
with  a  tunable  variation  of  the  Kraus  and  Turner  (1967)  one-dimensional  mixed 
layer  model  to  investigate  the  effect  of  mixed  layer  dynamics  on  the  thickness 
of  sea  ice  and  the  variation  in  the  ocean  heat  flux  into  the  ice.    The  tunable 
mixing  coefficients  in  the  mixed  layer  model  are  determined  by  solving  the 
Naval  Postgraduate  School  bulk  turbulence  closure  model  (Garwood,  1977). 
The  simplicity  of  the  Kraus-Turner  model  simplifies  the  investigation  of  the 
thermodynamic  and  dynamic  interactions  of  the  ocean  mixed  layer  and  sea  ice 
under  diverse  conditions.   Cases  that  need  to  be  studied  include  ice  response 
to  storms  under  varying  initial  temperature  and  salinity  profiles,  the  seasonal 
cycle  of  refreezing  and  melting,  and  longer-term  climatic  response  to  differing 
initial  thermodynamic  conditions  beneath  the  ice. 

Chapter  II  discusses  processes  involved  in  sea  ice  growth  and  ocean 
mixed  layer  dynamics.   In  Chapter  III  the  sea  ice  model  is  developed  and 
coupled  with  the  ocean  mixed  layer  model.   Chapter  IV  will  give  case  studies 
using  varying  initial  conditions  to  demonstrate  the  sensitivity  of  the  model  to 
different  forcing  and  the  interaction  between  dynamic  and  thermodynamic 


processes.   Chapter  V  summarizes  the  results  of  the  research,  and  makes 
conclusions  and  recommendations  for  further  study. 


II     THEORY 
A.      FORMATION  AND  GROWTH  OF  SEA  ICE 

The  presence  of  salt  in  seawater  depresses  both  the  freezing  temperature 
of  the  water  and  the  temperature  of  maximum  density,  compared  to  fresh  water. 
Usually,  the  salinity  of  seawater  is  greater  than  24.7  psu,  so  the  freezing 
temperature  is  higher  than  the  temperature  of  maximum  density.   Hence,  as  the 
surface  of  the  ocean  cools,  the  surface  water  density  increases  and  corrective 
mixing  results.  This  cooling  will  continue  until  the  mixed  layer  is  at  the  freezing 
point  (Weeks  and  Ackley,  1 982).   However,  the  underlying  thermodynamic 
structure  of  the  ocean  is  a  significant  factor  in  determining  the  possible  onset  of 
freezing  (Claes,  1990). 

When  the  mixed  layer  temperature  has  reached  the  freezing  point  any 
additional  heat  loss  causes  the  water  to  be  supercooled,  and  ice  begins  to 
form.   The  amount  of  supercooling  required  to  cause  ice  formation  is  on  the 
order  of  only  a  few  hundredths  of  a  degree  Celsius  and  leads  to  the  production 
of  Frazil  ice  (Maykut,  1985).   Depending  on  wind  and  wave  action  the  continued 
production  of  frazil  ice  results  in  "sheet  ice,H  grease  ice,  or  pancake  ice  (Weeks 
and  Ackley,  1982). 


Ice  formation  from  supercooling  will  continue  until  the  ice  forms  a 
continuous  sheet  that  tends  to  insulate  the  seawater  from  the  atmosphere.  As 
freezing  continues,  the  temperature  of  the   ice  surface  cools  from  near  the 
freezing  point  to  a  temperature  that  approaches  the  atmospheric  temperature. 
This  results  in  a  temperature  gradient  across  the  ice  and  a  conductive  heat  flux 
upward  through  the  ice,  removing  latent  heat  from  the  ocean.    Additional  ice 
growth  occurring  on  the  underside  of  the  sheet  (Maykut,  1985)  is  known  as 
congelation  growth  (Gow  and  Tucker,  1991).  This  new  ice  continues  to  develop 
and  may  form  different  types  of  ice  depending  on  environmental  conditions. 
These  ice  characteristics  are  important  for  understanding  and  predicting  ice 
growth,  and  ice  may  be  classified  by  stage  of  development  as  shown  in  Table 
2.1  (Naval  Polar  Oceanography  Center,  1984). 

The  amount  of  salt  initially  trapped  as  brine  in  sea  ice  is  a  function  of  the 
freezing  rate  of  the  sea  water  and  its  salinity.  The  makeup  and  distribution  of 
salt  in  the  ice  with  time  will  be  determined  by  the  temperature  of  the  ice  and  the 
phase  relationship  between  the  ice,  salt  and  brine  (Gow  and  Tucker,  1991). 
Figure  2.1    shows  idealized  vertical  salinity  profiles  of  sea  ice. 

Gow  and  Tucker  (1991)  summarized  the  important  features  of  the  salinity 
profiles  as  follows: 


High  salinities  at  the  top  and  bottom  of  thinner  ice  sheets,  leading  to 
strongly  C-shaped  salinity  profiles. 
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Figure  2.1     Idealized  Ice  Salinity  Profiles  in  Arctic  Sea  Ice  for  Various  Ice 
Thicknesses.   Curves  a  -  d  Represent  First  Year  Ice.   The  Remaining  Curves 
Illustrate  Profiles  Found  in  Multiyear  Ice  (Maykut,  1985) 


General  weakening  of  the  C-shaped  salinity  profile  accompanied  by 
progressive  decrease  in  bulk  salinity  with  increasing  thickness  and  age  of 
the  ice  sheet,  indicative  of  downward  migration  of  brine  through  the  ice 
and  its  ejection  into  the  underlying  ocean. 


Virtual  desalinization  of  the  upper  levels  of  First  year  ice  during  and  after 
its  transformation  to  multiyear  ice. 


The  shape  of  the  profiles  in  Figure  2.1  can  be  explained  by  four 
mechanisms  that  remove  brine  from  sea  ice:  migration  of  brine  pockets  through 
ice  crystals,  brine  expulsion,  brine  drainage,  and  flushing  (Untersteiner,  1968). 
Brine  drainage  is  probably  the  primary  mechanism  for  the  salinity  changes  in 
first  year  ice  (Maykut,  1985).   Cox  and  Weeks  (1974)  gathered  a  wide  variety  of 
data  to  show  the  dependence  of  average  salinity  upon  ice  thickness  in  the 
arctic.   Figure  2.2  shows  that  for  cold  (growing)  ice  there  is  a  steep,  linear 
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Figure  2„2  Average  Salinity  of  Sea  Ice  as  a  Function  of  Ice  Thickness  for  Cold 
Sea  Ice  Sampled  During  the  Growth  Season  (Cox  and  Weeks,  1974). 
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decrease  in  the  average  salinity  of  the  ice  up  to  an  ice  thickness  of  about  40 
cm,  and  then  a  sharp  decrease  in  the  slope  of  salinity  change  per  ice  thickness 
change.  The  abrupt  change  in  the  slope  is  most  likely  caused  by  a  shift  in  the 
desaiinization  mechanism  from  brine  expulsion  to  brine  drainage  (Gow  and 
Tucker,  1991).   The  high  initial  salinity  of  sea  ice  results  from  the  high  initial 


TABLE  2.1    STAGES  OF  ICE  DEVELOPMENT 


STAGE  OF  DEVELOPMENT 

THICKNESS 

New,  Frazil,  Slush,  Shuga,  Grease 

0  -  10  cm   (0-4  in) 

Dark  Nilas 

0  -  5  cm   (0-2  in) 

Light  Nilas 

5  -  10  cm  (2-4  in) 

Gray 

10-  15  cm   (4-6  in) 

Gray-White 

15 -30  cm   (6-12  in) 

|  First  Year  Thin 

30 -70  cm   (12-28  in) 

First  Year  Medium 

70-120  cm   (28-47  in) 

First  Year  Thick 

120 -200  cm   (47-80  in) 

Second  Year 

>  2  m   (>  80  in) 

Multiyear 

>  2  m  (>  80  in) 

growth  rate  of  the  ice  that  entraps  large  amounts  of  brine.  As  the  ice  thickens, 
the  growth  rate  of  the  ice  slows  and  less  brine  is  trapped.   The  effects  of  brine 
drainage  become  evident  and  result  in  a  decrease  in  the  salinity  of  the  ice 
(Maykut,1985). 
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For  simplicity  in  this  research,  the  ice  salinity  will  be  assumed  constant  in 
time,  but  the  Cox  and  Weeks  data  can  be  used  to  express  ice  salinity  as  a 
function  of  thickness. 

B.      MIXED  LAYER  DYNAMICS 

The  oceanic  planetary  boundary  layer  or  mixed  layer  is  the  fully  turbulent 
region  above  the  seasonal  thermocline.   Shear  stress  and  upward  buoyancy 
flux  provide  the  mechanical  energy  to  mix  salinity,  momentum  and  heat  in  the 
mixed  layer  (Garwood,  1977).   The  shear  stress  is  imparted  by  the  wind  to  the 
water  and  transmitted  downwards.  The  buoyancy  flux  results  from  either 
heating  from  below,  cooling  from  above  or  downward  surface  salinity  flux  which 
leads  to  instabilities  in  the  mixed  layer. 

The  mixed  layer  is  generally  considered  to  be  vertically  homogeneous  in 
temperature,  salinity  and  density  despite  the  presence  of  the  fine  scale 
microstructure    (Zilitinkevich  et  al.,  1979).  The  quasi-homogeneous  nature  of 
the  mixed  layer  has  led  to  it  being  described  as  a  "slab."    By  assuming  the 
mixed  layer  is  homogeneous  the  momentum  and  buoyancy  fluxes  only  need  to 
be  specified  at  the  surface  and  calculated  at  the  base  of  the  mixed  layer. 
However,  the  small  variations  in  the  mean  flow,  temperature  or  salinity  can  lead 
to  large  changes  in  the  turbulent  fluxes  (Garwood,  1977).   Zilitinkevich  et  al 
(1979)  show  that  the  vertical  heat  flux  from  a  0.01°  C  temperature  difference  in 
the  upper  1 0  meters  of  the  ocean  is  equivalent  to  the  heat  flux  generated  by  a 
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1°C  difference  between  the  sea  surface  and  the  atmosphere  at  a  height  of  10 
meters.     Figure  2.3  shows  a  temperature  and  salinity  profile  taken  at  an 
AID  J  EX  camp  (Bauer  et  al,  1980).   The  profiles  clearly  show  the  finescale 
structure  present  in  the  mixed  layer  that  is  not  shown  in  the  typical  temperature 
and  salinity  profiles  in  Figure  2.4.   Garwood  (1977)  showed  that  by  using  a 
vertically  integrated  (bulk)  model  of  the  mixed  layer  the  important  effects  of  the 
fine  scale  structure  of  the  mixed  layer  can  be  accounted  for  by  direct  estimation 
of  averaged  mixed  layer  turbulent  Kinetic  energy. 

The  Arctic  ocean  mixed  layer  has  some  unique  characteristics  due  to  its 
ice  cover.   Here  the  shear  stress  between  the  ice  and  ocean  provides  the 
mechanism  to  transfer  momentum  from  the  atmosphere  to  the  ocean. 
Typically,  the   ice-covered  mixed  layer  is  at  or  near  the  freezing  point  so  rt 
stores  almost  no  sensible  heat  (McPhee,  1986).   Because  the  density  of  water 
near  the  freezing  point  is  almost  entirely  a  function  of  salinity,  the  buoyancy  flux 
is  proportional  to  the  salt  flux  alone  (Morison  and  Smith,  1981). 

The  central  Arctic  is  divided  into  Canadian  and  Eurasian  basins  by  the 
Lomonosov  Ridge.   Coachman  and  Aagaard  (1974)  defined  the  three  water 
masses  in  the  central  Arctic:  Arctic  water  (0-200  m),  Atlantic  water  (200-900  m) 
and  bottom  water  (900  m-bottom).   The  Arctic  water  can  be  further  divided  into 
Arctic  surface  and  Arctic  subsurface  layer  waters.   Figure  2.4  shows  typical 
temperature  and  salinity  profiles  for  both  basins  of  the  arctic. 
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Caribou  Station  1  Oct  75  73.21N   139.06W 
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Figure  2.3  Temperature  and  Salinity  Profiles  Taken  With  a  CTD  During 
AIDJEX  from  Ice  Camp  Caribou  Showing  the  Finescale  Structure  That  is 
Typically  Present  in  the  Mixed  Layer  (Bauer  et  al,  1980). 


The  Arctic  surface  water  extends  from  the  surface  to  a  30  m-50  m  depth 
and  has  fairly  uniform  distributions  of  temperature  and  salinity  across  the  Arctic. 
The  variations  in  salinity  are  mainly  due  to  freshwater  runoff  from  land   and 
inflow  of  Atlantic  and  Pacific  Ocean  waters,  but  this  is  along  the  periphery  of 
the  basins.    The  temperature  is  at  or  near  the  freezing  point  throughout  the 
basins  except  for  areas  that  are   ice-free  for  extended  periods. 
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Figure  2.4  Vertical  Temperature  and  Salinity  Distribution  in  the  upper  400  m  as 
Observed  During  Winter  at  the  Three  Locations  Indicated  (after  Coachman  and 
Aagard,  1974). 

The  characteristics  of  the  Arctic  subsurface  layers  extend  from  the  surface 
waters  to  200  m  and  vary  between  the  two  basins.   In  the  Eurasian  basin  the 
salinity  increases  with  depth  down  to  about  200  m  where  it  reaches  a  value  of 
34.9  to  35.0  psu,  and  below  there  it  is  nearly  isohaline.   The  temperature 
remains  isothermal 
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down  to  about  1 00  m  and  then  increases  to  above  0°  C  due  to  the  presence  of 
the  Atlantic  water.   In  the  Canadian  basin,  the  halocline  is  similar  to  the 
Eurasian  basins  except  there  is  not  as  rapid  an  increase  in  salinity  with  depth. 
The  major  Canadian  basin  difference  is  in  the  temperature  profile  due  to  the 
presence  of  Bering  sea  water,  which  results  in  a  temperature  maximum  of 
approximately  -1°  C  at  75  to  100  m  depth.  The  temperature  then  decreases  to 
-1.5°  C  at  about  150  m  and  remains  nearly  isothermal  down  to  the  Atlantic 
water.  The  presence  of  the  Bering  sea  water  is  critical  in  determining  the 
oceanic  heat  flux  into  the  ice  because  it  separates  the  surface  layer  from  the 
warm  (>0°C)  Atlantic  water,  as  shown  by  Morison  and  Smith  (1981)  from 
observations  on  ice  island  T-3. 
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III.     MODEL  DEVELOPMENT 

A.      THE  ICE  GROWTH  MODEL 

The  physical  and  dynamical  properties  of  sea  ice  have  been  studied  for 
many  years  but  are  still  not  completely  understood.  However,  the  most 
important  aspects  of  sea  ice  have  been  extensively  investigated  and 
parameterizations  and  numerical  models  of  varying  degrees  of  sophistication 
have  been  developed  to  explain  the  behavior  of  sea  ice.   For  this  model  only 
the  most  important  physical  processes  will  be  included  and  will  be  treated  in  as 
simplified  a  manner  as  possible.   The  simplifications  are  not  ad  hoc  but  are 
based  on  the  nature  of  the  physical  process  involved  and  the  extent  it  is 
understood.   This  method  may  lead  to  minor  quantitative  differences  compared 
to  more  sophisticated  models,  but  the  goal  here  is  to  emphasize  the  role  of 
mixed  layer  dynamics  in  limiting  sea-ice  growth.   Hence,  complexity  attributed  to 
other  processes  is  minimized. 

1.      Longwave  Radiation  Fluxes 

The  first  step  is  to  linearize  the  Stefan-Boltzman  Law  near  the 
freezing  point  to  determine  the  upward  radiant  flux  leaving  the   ice.   Following 
Thorndike  (1992)  the  Stefan-Boltzman  Law  can  be  written  as: 
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F*r°T*r°Vf>Tj  (3.D 


up 

where: 

a  =  5.7  x  10"8  Wm"2K^   is  the  Stefan  Boltzman  constant 
T,=273  K  is  the  freezing  point  of  fresh  water 
T,^  is  the  surface  temperature  in  °C 
Fup  is  the  radient  exitance  in  Wm"2 


Linearizing  equation  3.1  leads  to: 


F^A+BT*.  <3.2) 


where 

A  =  o(Tf)4  =  320  Wm"2 
B  =  4o(Tf)3  =  4.6  Wm"2  K"1 


Using  a  surface  temperature  of  -40°C  with  equation  3.1  gives  a  radiation 
flux  of  168  Wm'2  while  equation  3.2  gives  a  radiation  flux  of  136  Wm'2  which  is 
a  difference  of  19%.  The  downward  long  wave  radiation  is  paramaterized  ( 
Maykut,  1986)  as: 

F^e-oTl  (3.3) 

where: 

e*  is  the  effective  ernissivity  of  the  atmosphere 
o  is  the  Stefan  Boltzman  constant 

Tall  is  the  air  temperature  at  2  m  and  will  be  specified  to  simulate  various 
atmospheric  conditions 

Maykut  and  Church  (1973)  developed  an  expression  for  e*  based  on 

observations  at  Barrow,  Alaska  which  is: 
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€•=0.7855(1  +0.2232C2*75) 

where: 

C  is  the  fractional  cloud  cover 


2.      Sensible  Heat  flux 

Sensible  heat  flux  in  the  atmospheric  surface  layer  above  the  ice  is 
the  result  of  turbulent  heat  conduction  between  the  atmosphere  and  the  ice 
surface.  Although  the  sensible  heat  flux  is  often  ignored  in  climate  studies  of 
the  central  Arctic,   it  can  be  very  important  on  shorter  time  scales.   Figure  3.1 
shows  the  variations  of  the  sensible  heat  flux  with  ice  thickness  and  season.   In 
early  March  the  sensible  heat  flux  to  the  atmosphere  exceeds  350  Wm'2  for  5 
cm  thick  ice.  The  importance  of  the  sensible  heat  flux  decreases  with 
increasing  thickness  of  ice,  and  the  sensible  heat  flux  for  0.8  m  thick  ice  is 
about  40  Wm'2.   The  sensible  heat  flux  has  been  para  mate  rized  by  Maykut 
(1986)  as: 

where 

p8  is  the  density  of  air 

c^  is  the  specific  heat  of  air 

cs  is  the  bulk  transfer  coefficient  for  sensible  heat 

u10,v10  is  the  wind  speed  at  10m 
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TBy  is  the  air  temperature  at  2m 

T^  is  the  temperature  of  the  surface  of  the  ice 

Andreas  (1987)  has  found  the  value  for  Cs  to  vary  between  1.0  x  10 

3  and  1.5  x  10"3.  For  this  model  the  value  of  Cs  will  be  fixed  at  1.25  x  10"3.  The 

air  temperature  at  2  meters  will  be  specified  to  simulate  various  atmospheric 

conditions. 

3.      Conductive  Heat  Flux 

After  the  summer  melting  season  the  surface  of  the  ice  begins  to 
cool  while  the  bottom  of  the  ice  remains  near  the  freezing  point  of  the  adjacent 
seawater.  When  the  ice  surface  becomes  colder  than  the  ice  bottom  an 
upward  conduction  of  heat  results. 

This  conductive  heat  flux  is  the  amount  of  heat  reaching  the  surface 
and  depends  on  the  temperature  gradient  across  the  ice  and  the  thermal 
conductivity  of  the  ice.   Maykut  and  Untersteiner  (1971)  paramaterized  the 
conductive  heat  flux  as: 


F=-K?I  (3.5) 

c         dz 


Where: 


Fc  is  the  conductive  heat  flux 

K  is  the  thermal  conductivity  of  ice 

dt 

—  is  the  temperature  gradient  of  the  ice 
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Figure  3.1  Seasonal  Variation  in  the  Sensible  Heat  Flux  as  a  Function  of  Ice 
Thickness  (m)  (Maykut,  1986). 

The  thermal  conductivity  (K)  of  young  ice  can  vary  by  as  much  as 
50%  (Maykut,  1986).    However,  following  Thorndike  (1992)  K  is  considered  to 
be  constant  of  2  Wm1  K"\   The  temperature  profile  of  the  ice  will  be  considered 
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linear,  which  Maykut  (1978)  showed  to  be  a  good  approximation  for  young  ice 
and  allows  equation  3.5  to  be  written: 


FC=K  f     **  (3.6) 

c  H 


where. 


H  is  the  thickness  of  the  ice 

T^  is  the  temperature  of  the  top  of  the  ice 

Tf  is  the  temperature  of  the  bottom  of  the  ice  using  equation  3.7 

a.      Freezing  point  of  sea  water 

Neglecting  pressure  effects,  the  freezing  point  temperature  of 
sea  water  is  given  by  equation  (15)  (Millero  and  Leung,  1976): 


T^aoS^+aiS^  +  azS?  <3-7> 


where: 


a0  =  -0.0575 

a,  =  1.710523  x  10  3 

a2  =  -2.154996  x  10^ 


s,  is  the  salinity  of  the  mixed  layer 

4.      Ice  Motion 

The  equation  of  motion  for  sea  ice  balances  ice  acceleration  against 
air  and  water  stresses,  coriolis  force,  internal  ice  stress  and  the  pressure 
gradient  force.  Thorndike  and  Colony  (1982)  showed  that  over  70%  of  the 
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variance  of  the  ice  motion  can  be  explained  by  the  geostrophic  wind  and 
developed  a  free  drift  model  that  neglected  the  internal  ice  stresses.   Hibler  and 
Tucker  (1979)   show  that  in  the  short  term  the  effects  of  the  pressure  gradient 
force  is  small  compared  to  the  coriolis  force  and  wind  and  water  stresses  and  is 
therefore  neglected  in  this  model.   In  any  case,  the  geostrophic  flow  may  be 
linearly  separated  from  the  force  balance,  and  the  equations  for  unsteady  ice 
motion  become: 

ZZLmfr*    *""      x"  (3.8) 


dt     '     P/  H 


dv,  tv    +*v 

—±=-fU,+    ¥m     Yw  (3*) 

dt  '        9,H 


where: 


u„  v,  are  the  x  and  y  components  of  ice  velocity 

f  is  the   Coriolis  Force 

iXBtm,  xy8tm   is  the  stress  imparted  on  the  ice  from  the  wind. 

xxw,  xyw  is  the  stress  imparted  on  the  ice  from  the  ocean. 

p,  is  the  density  of  the  ice. 

H  is  the  thickness  of  the  ice. 

a.      Wind  and  Ocean  stresses 

The  wind  stress  on  the  ice  can  be  paramaterized  as: 


ta^VJ^riWrV  (310) 
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V«,=P.^.(Vio-^Mo-^z+(^io-^2  (311) 

While  the  ocean  stress  on  the  ice  can  be  parameterized  as: 

V  t>.C„lUrUjJ(UrUj/'+(VrVjr  (3«> 

V."  PwCDl{Vrvy{UrUjfHVrVj  (3-13) 


where: 


pa  is  the  density  of  air 

pw  is  the  density  of  sea  water 

C^  is  the  ice  surface  drag  coefficient 

CDI  is  the  ice  bottom  drag  coefficient 

u10,  v10  are  the  x  and  y  components  of  the  wind  speed  at  10m 

u„  v,  are  the  x  and  y  components  of  the  ice  velocity 

u*.  vw  are  the  x  and  y  components  of  the  ocean  mixed  layer  velocity 


Guest  and  Davidson  (1991)  measured  the  surface  drag 
coefficient  (  C^  )  of  various  ice  types.   For  this  model  a  constant  value  of  2.3  x 
1 0 3  will  be  used  for  C^  which  corresponds  to  Guest  and  Davidson's  median 
value  for  smooth  young  ice.   McPhee  (1990)  emphasized  the  difficulty  in 
determining  a  bottom  drag  coefficient  and  the  wide  range  of  values  of  C,^  that 
have  been  obtained.   However,  this  model  uses  a  constant  value  of  5.4  x  10"3 
for  the  bottom  drag  coefficient  following  McPhee  (1979). 
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5.  Oceanic  Heat  Flux 

The  ocean  heat  flux  is  the  result  of  wanner  water  being  entrained 
into  the  mixed  layer  due  to  turbulent  mixing  and  is  paramaterized  by  Chu  and 
Garwood  (1988)  as: 

Fws'PwOPwW9LT  (3.14) 

where: 

Fw  is  the  ocean  heat  flux  due  to  entrainment 
pw  is  the  density  of  sea  water 
we  is  the  entrainment  rate 

AT  is  the  temperature  difference  between  the  ocean  mixed  layer  and  the 
water  below. 

The  entrainment  rate  we  is  given  in  equation  3.28  and  is  determined  by  the  one- 
dimensional  mixed  layer  model  which  is  discussed  in  section  B  of  Chapter  III. 

6.  Ice  Growth 

The  rate  of  ice  growth  or  melting  will  be  determined  by  the  balance 
of  fluxes  into  and  out  of  the  ice.     Latent  heat  is  released  at  the  bottom  of  the 
ice  by  changing  water  into  ice  when  the  conductive  heat  flux  upward  into  the  ice 
exceeds  the  ocean  heat  flux  upward  into  the  mixed  layer.  The  ice  will  melt  If 
the  ocean  heat  flux  is  greater  than  the  conductive  heat  flux  (Thorndike,  1 992). 
The  latent  heat  flux  due  to  the  change  in  ice  thickness  is  given  as: 
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'     at 


Where: 

F,  is  the  latent  heat  flux  due  to  change  in  ice  thickness 

L  is  the  latent  heat  of  fusion  for  ice 

dH/dt  is  the  rate  of  change  in  ice  thickness 

The  latent  heat  of  fusion  is  considered  to  be  constant  at  3  x  108  J  m"3 

The  balance  of  heat  fluxes  is  represented  in  Figure  3.2.   Summing 

the  heat  fluxes  at  the  surface  of  the  ice  yields: 


Fl*+Fc+F.=Fv  <3-15> 


and  summing  the  heat  fluxes  at  the  bottom  of  the  ice  yields: 


iM  =  Fc-Fw  (3.16) 

at 


Using  equations  3.2,  3.4  and  3.6  with  equation  3.15  gives: 


Fu»+K  ^^-P^W  C7"*- TJ  +A+Br*    (317) 


Using  equations  3.6  and  3.14  with  equation  3.16  gives: 
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^'K®T^+'M*T 


(3.18) 


Solving  equation  3.17  for  the  surface  temperature  yields: 


FLW-A+pacpc8Ju?0+v?0Ta 
'sic* ;  .     .    j—  + 


KT, 


13.19) 


*+P^P/W^Wo+-£       BH+  PacpCjJu?o+v?0  H+K 


Figure  3.2    The  Heat  Fluxes  That  Determine  The  Ice  Growth  and  Temperature 
( from  Thorndike,  1992) 
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Finally,  solving  equation  3.18  for  ice  growth  gives: 

dH^K(TrTsfc)  +  l_{  A  J)  (3.20) 

B.      OCEAN  MIXED  LAYER  MODEL 

The  mixed  layer  model  used  here  is  a  tunable,  one  dimensional  Kraus  and 
Turner  (1967)  type  mixed  layer  model.  The  model  assumptions  following  Niiler 
and  Kraus  (1977)  are: 


1 .  The  mean  temperature,  salinity  and  horizontal  velocity  are  assumed  to 
be  uniform  within  the  mixed  layer. 

2.  On  the  depth  and  time  scales  of  the  model,  temperature,  salinity  and 
horizontal  velocity  discontinuities  can  exist  across  the  lower  boundary  of 
the  mixed  layer. 

3.  Temperature  changes  associated  with  frictional  dissipation  are  neglected. 

4.  Wind   and  ice  stresses,  turbulent  kinetic  energy  and  viscous  dissipation 
are  considered  to  be  horizontally  uniform. 


Additional  assumptions  are: 


5.  The  temperature  of  the  mixed  layer  remains  at  the  freezing  point,  thus 
not  allowing  storage  of  heat  by  the  mixed  layer  and  neglecting  the  effect 
of  supercooling  at  the  water  surface. 

6.  Loss  of  turbulence  to  internal  waves  is  ignored. 

7.  Only  the  case  of  deepening  mixed  layer  will  be  considered  at  the  present 
time. 
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1.      Turbulent  Kinetic  Energy  Equation 

The  model  is  derived  from  the  turbulent  kinetic  energy  (TKE) 
equation,  equation  3.21  (Garwood,  1977): 


1  a(i//2+^5+n^/2)_ 

dz         dz. 

-  w'*'g- a 

Po       dz 

w' 

u^+v^+w'*   p'\ 

2            dt 

1         2         *  Pol 

c 

(3.21) 

i 

II 

III 

IV 

V 

From  Stull  (1988): 


•  Term  I   is  the  local  storage  tendency  of  TKE. 

•  Term  II   is  the  shear  production  of  TKE  due  to  vertical  shear  of  the 
horizontal  mean  flow. 

•  Term  III  is  the  buoyant  production  or  consumption  of  TKE  depending  on 
whether  the  buoyancy  flux  is  positive  (upwards)  or  negative  (downward). 

•  Term  IV  is  turbulent  transport  of  TKE  by  turbulent  eddies  and  the  pressure 
transport  of  TKE  due  to  pressure  perturbations  and  can  be  positive  or 
negative. 

•  Term  V  is  the  viscous  dissipation  of  TKE  due  to  the  conversion  of  TKE 
into  heat  and  will  always  be  negative. 


The  TKE  budget  of  the  ocean  can  be  assumed  to  be  in  steady  state 
because  of  the  very  short  time  scale  for  dissipation.   This  allows  Term  I  of 
equation  3.21  to  be  assumed  equal  to  zero. 
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2.      Vertically  Integrated  TKE  Equation 

Equation  3.21  is  now  vertically  integrated  over  the  depth  of  the  mixed 
layer  and,  following  Claes  (1990)  ,  gives  the  bulk  form  of  the  TKE  equation: 


0 =-^-  +[ag  7V(0)  -  $g  5  V(0)]  +[aflf  T'w'(-h)  -  fig  sV(-jf?jf 2) 
h 


ii  hi 


where: 

c1  is  a  constant  of  proportionality 

u.  is  the  water  surface  friction  velocity 

h  is  the  depth  of  the  mixed  layer 

a  is  the  thermal  expansion  coefficient 

(3  is  the  salinity  contraction  coefficient 

g  is  gravity 

7"  V(0)  is  the  temperature  flux  at  the  top  of  the  mixed  layer 

sV(0)    is  the  salinity  flux  at  the  top  of  the  mixed  layer 

T'w'{-fy  is  the  temperature  flux  at  the  base  of  the  mixed  layer 

s'w'(-h)  is  the  salinity  flux  at  the  bottom  of  the  mixed  layer 

Because  the  temperature  of  the  mixed  layer  is  kept  at  the  freezing  point,  any 
addition  or  loss  of  heat  will  go  to  melting  or  freezing  ice.   Hence  the  water 


surface  temperature  flux  agT'w'(Q)  is  equal  to  zero 
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a.      Ice  stress  and  TKE  transport 

Term  I  in  equation  3.22  represents  the  shear  production  minus 
viscous  dissipation  integrated  over  the  full  depth  of  the  mixed  layer.   The 
constant  c,  is  reduced  exponentially  on  a  vertical  scale  of  u./f  (Chu  and 
Garwood,  1988)  and  is  approximately: 


c,«2* 


i-fh/u.) 


where: 


f  is  the  coriolis   parameter 

h  is  the  depth  of  the  mixed  layer 

u.  is  the  water  surface  friction  velocity. 

c1   can  be  tuned  to  match  observations.   The  friction  velocity  (u. )  is  represented 

by: 


2  2 

xi  71 


T>/2 


where: 

xxl ,  Ty,  are  the  stresses  imparted  on  the  water  from  the  ice  . 
pw  is  the  density  of  the  sea  water. 

The  ice  stresses  x,,,  and  Xy,  are  equal  but  opposite  in  sign  to  the  stress  imparted 

on  the  ice  from  the  water,  given  in  equations  3.12  and  3.13. 

b.       Buoyancy  fluxes 

The  buoyancy  flux  at  the  surface  of  the  mixed  layer  will  be 
determined  by  the  salinity  flux  at  the  top  of  the  mixed  layer  due  to  ice  formation 
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or  melting,  while  the  buoyancy  flux  at  the  base  of  the  mixed  layer  will  be 
determined  by  the  heat  and  salinity  fluxes  due  to  entrapment  of  water  from 
below  the  mixed  layer. 


(1)     Surface  buoyancy  flux.   The  s  V(0)  term  of  equation  3.22 
is  the  surface  salinity  flux  that  results  from  the  freezing  or  melting  of  ice  and  is: 

5V(0)  =  -^-(Sr  S,)4r  (3.23) 

9W  dt 

where: 

p,  is  the  density  of  the  ice 

pw  is  the  density  of  the  sea  water 

s1  is  the  salinity  of  the  mixed  layer 

s,  is  the  salinity  of  the  ice 

dH/dt  is  the  growth  rate  of  the  ice 


Term  II  of  equation  3.22  is  the  surface  buoyancy  flux  b'w'(0)  and  using 
equation  3.23  can  be  written  as: 

P//-    ^dH 


b'w'W^gSLft-Si)^  (3.24) 

9w  dt 


(2)     Buoyancy  flux  at  the  base  of  the  mixed  layer.  The 


T  w\-h)  and  s  w'(-h)    terms  of  equation  3.22  are  the  heat  and  salinity  fluxes 

at  the  base  of  the  mixed  layer  that  result  from  the  entrainment  of  water,  with 
different  temperature  and  salinity,  from  below  the  mixed  layer .   Garwood 
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(1977)  develops  the  jump  condition  for  the  turbulent  fluxes  at  the  base  of  the 
mixed  layer  by  integrating  across  the  entrapment  zone  and  assuming  a 
negligible  amount  of  heat  stored  in  the  entrainment  zone.   This  allows  the 
temperature  flux  at  the  base  of  the  mixed  layer  to  be  written  as: 


T'w'(-h)  =  -LTwe  <3-25> 

where: 

A  7=  Tr  T2 


AT  is  the  temperature  jump  between  the  mixed  layer  temperature  and  the 
deep  water  temperature 

Tf  is  the  mixed  layer  temperature  which  is  at  the  freezing  point 
T2  is  the  deep  water  temperature 


and  the  salinity  flux  to  be  written  as: 


$/w,{-h)^-Asw6  <3-26) 


where: 


A5  =  51-^ 
where  : 


As  is  the  salinity  jump  between  the  mixed  layer  salinity  and  the  deep  water 
salinity. 

s,  is  the  mixed  layer  salinity 
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s2  is  the  deep  ocean  salinity 
we  is  the  entrain  ment  velocity 


Term  III  of  equation  3.22  is  the  buoyancy  flux  b'w'(-h)  at  the  base 
of  the  mixed  layer.   Using  equations  3.25  and  3.26  can  be  written  as: 


b'w'i-Q-agrw'i-fy  -$gs'w'(-h) 


(3.27) 


with  : 


Ab=ag&T-pgAs 


Ab  is  the  buoyancy  jump  between  the  mixed  layer  and  the  deep  water 

3.       Entralnment  Velocity 

The  entrainment  velocity  can  be  obtained  from  equation  3.22 
combined  with  equations  3.24  and  3.27  and  an  added  tuning  constant  c2 ,  which 
is  a  function  of  stability,  to  yield: 


h  pw  at 


(3.28) 


Solving  equation  3.28  for  the  entrainment  velocity,  We  gives: 


w,~ 


Ab 


h  pw  at 


or 


(3.29) 
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w^ 


6    ugbT-pgAs 


h  pw  dt 


(3.29) 


4.      Mixed  Layer  Motion 

The  equations  of  motion  for  the  mixed  layer  are  determined  by 
integrating  the  mean  momentum  equation  across  the  depth  of  the  mixed  layer 
(Garwood,  1977),  resulting  in  equations  3.30  and  3.31 


duw     *  xx>       UWWB  „_„ 


ar      •"  p^/7     /? 


^.-ftf   +JjLJ5^-flL|/  (3.31) 

dt  w    pwh      h       3  w 


where: 

uw ,  vw  are  the  horizontal  velocities  of  the  mixed  layer 

f  is  the  coriolis  force 

h  is  the  depth  of  the  mixed  layer 

pw  is  the  density  of  the  mixed  layer 

we  is  the  entrainment  rate 

x^ ,  Ty,  are  the  stresses  imparted  on  the  mixed  layer  by  the  ice  motion 

c3  is  a  linear  damping  constant. 

The  linear  damping  term  (Cg)  has  been  added  to  equations  3.30  and 

3.31  to  model  the  momentum  dispersion  effect  (Pollard  and  Millard,  1970)  and 

can  be  tuned  to  match  observations.   Typically,  Cg  *  (2  days)"1  . 
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5.      Mixed  Layer  Salinity 

The  salinity  of  the  mixed  layer  will  vary  as  a  result  of  the  salinity 
fluxes  given  by  equations  3.2  and  3.26.  The  salinity  equation  is  determined  by 
integrating  the  mean  salinity  equation  across  the  depth  of  the  mixed  layer  and 
is  given  as: 


dt     h 


Pw  w 


(3.32) 
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IV.     MODEL  RESULTS 


A.      MODEL  INITIAL  CONDITIONS 

The  closely  coupled  nature  of  ice  growth  with  atmospheric  and  oceanic 
forcings  does  not  lend  itself  to  a  sensitivity  study  where  only  one  parameter  is 
varied  while  all  others  are  held  constant.    However,  by  varying  initial  conditions 
and  boundary  conditions  in  the  model,  the  sensitivity  of  the  model  to  the  various 
parameters  can  be  determined.  The  model  boundary  conditions  are  the  air 
temperature,  wind  speed,  and  cloudiness.  The  initial  conditions  are  ice  velocity, 
mixed  layer  velocity,  initial  mixed  layer  salinity,  the  deep  layer  temperature,  the 
deep  layer  salinity  and  the  initial  mixed  layer  depth.  Table  4.1  lists  the  case 
solutions  and  the  parameter  that  was  varied.   Output  for  all  model  cases  is 
shown  in  the  appendix.  The  first  case  is  the  standard  against  which  the  effect 
of  the  variations  will  be  considered  and  is  based  on  the  temperature  and 
salinity  profiles  shown  in  Figure  2.1.      For  all  runs  the  initial  velocity  of  the  ice 
and  mixed  layer  are  kept  constant  as  is  the  cloudiness.   One  of  the  major  goals 
of  this  research  is  to  study  the  variation  in  the  oceanic  heat  flux  and  its  effect 
on  the  growth  rate  of  ice. 
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1.  VARYING  WIND  SPEED 

Calm  winds  were  not  considered  in  this  study,  although  this  may  be 
an  important  situation.   When  the  winds  calm,  the  mixed  layer  will  shallow, 
entrapment  will  cease,  and  there  will  be  no  ocean  heat  flux  into  the  ice. 

The  standard  case  wind  speed  is  chosen  to  be  approximately  7ms 1 
while  in  case  2  the  wind  speed  is  doubled  to  14ms'1  .  The  increased  wind 
speed  has  two  effects  on  the  ice  growth  rate.   First,  the  sensible  heat  flux  is 
increased  which  decreases  the  surface  temperature  of  the  ice.   Second,  as  can 
be  seen  from  equation  3.29,  the  mechanical  mixing  (u.)  term  will  increase  eight- 
fold when  the  wind  speed  (and  u.)  is  increased  by  a  factor  of  two.  This  leads  to 
a  large  entrapment  rate  (we )  and  a  heat  flux  that  exceeds  the  conductive  heat 
flux  (Fc )  in  less  than  four  days  and  begins  to  melt  the  ice. 

2.  VARYING  AIR  TEMPERATURE 

When  the  air  temperature  is  increased  from  -30.0°  C  to  -10.0°  C,  the 
incoming  net  longwave  radiation  is  increased  to  220  Wm 2.  This  increased  heat 
into  the  ice  increases  the  surface  temperature  of  the  ice  to  the  point  that  the 
sensible  heat  flux  becomes  positive  which  also  leads  to  a  warmer  surface 
temperature.   The  warmer  surface  temperature  decreases  the  conductive  heat 
flux  and  slows  the  growth  rate  of  the  ice.   The  slower  ice  growth  has  two  effects 
on  the  ocean  heat  flux.   First,  the  surface  buoyancy  flux  is  decreased  which 
slows  the  entrapment  rate  (we )  and  second,  there  is  less  salt  being  ejected  into 
the  mixed  layer.   Since  the  freezing  temperature  is  determined  by  the  mixed 
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layer  salinity,  the  temperature  jump  (AT)  is  also  decreased.   This  results  in  a 
lower  ocean  heat  flux,  however,  the  ocean  heat  flux  is  actually  a  larger 
percentage  of  the  conductive  heat  flux  in  this  case. 

3.  VARYING  THE  ICE  THICKNESS 

Decreasing  the  initial  ice  thickness  (H)  to  0.05  m  causes  the  surface 
temperature  to  increase  but  still  results  in  almost  doubling  the  conductive  heat 
flux    (Fc ).   This  large  conductive  heat  flux  rapidly  decreases  as  the  ice 
thickens  because  the  temperature  gradient  across  the  ice  changes  much  more 
slowly  than  does  the  change  in  ice  thickness.   The  rapid  ice  growth  causes  an 
increase  in  the  entrainment  rate  and  the  ocean  heat  flux.   The  increase  in 
entrainment  rate  manifests  itself  in  the  increased  depth  of  the  mixed  layer  (h). 

Increasing  the  initial  ice  thickness  (H)  to  2m  allows  the  surface 
temperature  to  decrease  but  results  in  almost  a  50%  decrease  in  the 
conductive  heat  flux  and  a  slowing  in  the  ice  growth  rate.   This  slower  ice 
growth  rate  results  in  a  slower  change  in  the  buoyancy  jump  (Ab)  between  the 
mixed  layer  and  the  deep  layer  because  the  mixed  layer  is  now  fresher  than  in 
the  standard  case. 

4.  VARYING  THE  MIXED  LAYER  DEPTH 

Decreasing  the  initial  mixed  layer  depth  (h)  sharply  increases  the 
entrainment  rate  for  two  reasons.   First,  the  mechanical  mixing  is  inversely 
proportional  to  the  mixed  layer  depth  (h)  times  Ab.   Second,  the  surface 
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buoyancy  flux  is  increased  due  to  the  rapid  increase  in  mixed  layer  salinity 
which  is  also  inversely  proportional  to  the  mixed  layer  depth  as  shown  in 
equation  3.32.   Increasing  the  initial  mixed  layer  depth  to  45m  likewise 
decreases  the  ocean  heat  flux  as  the  entrainment  rate  (we)  becomes 
proportionally  smaller  to  the  entrainment  rate  in  the  standard  case. 

5.  VARYING  THE  DEEP  LAYER  SALINITY 

Increasing  the  deep  layer  salinity  (s2)  to  31.25  psu  reduces  the 
entrainment  velocity  (we )  because  the  deep  water  is  now  denser  and  requires 
more  work  to  lift  it  into  the  mixed  layer.  The  salinity  jump  As  does  decrease  at 
a  more  rapid  rate  due  to  the  increased  ice  growth,  but  it  is  not  enough  to  offset 
the  larger  initial  As. 

Decreasing  the  deep  layer  salinity  to  30.35  psu  results  in  a  large 
increase  in  the  entrainment  rate  (we )  because  of  the  decreased  buoyancy  jump 
(Ab).  The  decrease  in  the  ice  growth  and  consequent  decrease  in  the  surface 
buoyancy  flux  is  very  small  compared  to  the  increased  entrainment  rate. 

6.  VARYING  THE  DEEP  LAYER  TEMPERATURE 

Increasing  the  deep  layer  temperature  (T2)  to  -1.2°C  increases  the 
ocean  heat  flux  by  increasing  the  temperature  jump  (AT).   However,  this 
produces  only  a  slight  change  in  the  entrainment  velocity  because  the  thermal 
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expansion  coefficient  (a)  is  so  much  smaller  than  the  salinity  contraction 
coefficient  (p). 

The  effect  of  decreasing  T2  so  that  the  ocean  is  essentially 
isothermal  is  not  shown  here,  but  would  result  in  a  vanishing  ocean  heat  flux. 

7.       VARYING  THE  MIXED  LAYER  SALINITY. 

Increasing  the  initial  mixed  layer  salinity  (s1 )  to  30.35  psu  not  only 
decreases  the  salinity  jump  (As)  but  it  also  increases  the  temperature  jump  (AT) 
which  leads  to  larger  entrapment  rates  (we )  and  consequently  a  larger  ocean 
heat  flux  (Fw ).    A  decrease  in  salinity  leads  to  a  similar  but  opposite  result, 
although  this  is  not  shown  here. 
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TABLE  4.1   STANDARD  RUN  VALUES  AND  VARIED  PARAMETER  VALUES 


RUN# 

PARAMETER  VARIED 

PARAMETER  VALUE 

Run  1 

Standard  Case 

Ta=-30.0°C      u10=  5.0  ms1 

v10=  5.0  ms'1 

s1  =  29.9  psu    Sj=  30.8  psu 

T2=-1.4°C 

H  =  1.0  m   h=30.0m 

Run  2 

Wind  Speed 

u10=  10.0  ms"1   v10=  10.0  ms"1 

Run  3 

Air  Temperature 

Ta=-10.0°C 

Run  4 

ice  Thickness 

H  =  0.5  m 

Run  5 

Ice  Thickness 

H  =  2.0  m 

Run  6 

Mixed  Layer  Depth 

h  =  15.0  m 

Run  7 

Mixed  Layer  Depth 

h  =  45.0  m 

Run  8 

Salinity  Jump,  As 

s2  =  31 .25  psu 

Run  9 

Salinity  Jump,  As 

Sj  =  30.35  psu 

Run  10 

Temperature  Jump, 
AT 

T2=-1.2°C 

Run  11 

Salinity  Jump,  As 
Temperature  Jump, 
AT 

s1  =  30.35  psu 

40 


V.   CONCLUSIONS 

A.       SUMMARY 

This  study  has  investigated  the  effects  of  mixed  layer  dynamics  on  the 
growth  of  sea  ice  using  a  coupled  ice-mixed  layer  model.    The  ice  is  modeled 
as  undeformed  ice  with  a  linear  temperature  profile  and  constant  conductivity. 
The  ocean  is  modeled  as  a  one  dimensional,  two  layer  system.   Important 
prognostic  variables  include  ice  thickness,  mixed  layer  salinity,  and  mixed  layer 
depth.  As  long  as  ice  is  present,  the  mixed  layer  temperature  is  held  at  the 
freezing  point,  which  is  dependent  upon  the  mixed  layer  salinity.   The  surface 
temperature  of  the  ice  varies  in  response  to  the  net  surface  heat  fluxes. 

Model  solutions  revealed  the  response  of  the  model  to  initial  conditions, 
varying  one  parameter  at  a  time.  These  parameters  included  wind  stress,  air 
temperature,  mixed  layer  depth,  ice  thickness,  and  temperature  and  salinity 
jumps  between  the  mixed  layer  and  deep  layer.  The  model  simulations 
showed: 


•  There  is  a  negative  feed  back  between  the  ice  growth  rate  and  the  mixed 
layer  entrainment  rate. 

•  Increasing  wind  stress  leads  to  a  rapid  increase  in  the  entrainment  rate 
into  the  mixed  layer  and  an  upward  ocean  flux  that  exceeds  the 
conductive  heat  flux. 
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•  The  mixed  layer  entrainment  rate  is  strongly  influenced  by  the  ice  growth 
rate  through  its  effects  on  the  surface  buoyancy  flux,  mixed  layer 
temperature  and  salinity  and  subsequent  effect  on  the  buoyancy  jump  Ab. 

•  Varying  the  deep  layer  salinity  has  a  greater  effect  on  the  ocean  heat  flux 
than  does  varying  the  deep  layer  temperature. 

•  Increasing  the  mixed  layer  salinity  leads  to  a  large  increase  in  the  ocean 
heat  flux  through  the  combined  effects  of  increasing  the  temperature  jump 
AT  and  decreasing  the  buoyancy  jump  Ab. 

•  Increasing  the  mixed  layer  depth  decreases  the  entrainment  rate  by 
decreasing  the  effects  of  wind  stress  and  decreasing  the  surface 
buoyancy  flux.   Decreasing  the  mixed  layer  depth  increases  the 
entrainment  rate  for  the  same  reasons. 

•  The  large  range  of  heat  fluxes  precludes  modelling  the  heat  flux  as  a 
constant  or  using  more  simplified  physics  to  determine  the  ocean  heat 
flux. 

B.       RECOMMENDATIONS 

The  model  shows  the  importance  of  mixed  layer  dynamics  on  ice  growth 

during  the  freezing  cycle  and  under  active  mixed  layer  deepening.  This 

research  has  highlighted  several  areas  of  improvement  and  areas  of  further 

study: 


•  Incorporating  the  three  dimensional  Naval  Postgraduate  School  mixed 
layer  model  (Adamec  et  al.,  1981)  with  the  ice  growth  model  to  allow  for 
deepening  and  shallowing  of  the  mixed  layer  and  advection. 

•  The  effect  of  storms  and  varying  wind  stresses  on  the  ocean  heat  flux. 

•  Determination  of  the  equilibrium  thickness  of  sea  ice  through  an  entire 
year. 


42 


APPENDIX      MODEL  RUN  RESULTS 


Run  1: 


INITIAL  CONDITIONS 


u,    =0.1  ms"1 
s,  =9.0  psu 
v,    =0.1  ms'1 
i^  =  3.0x1 0"3  ms"1 
vw  =  3.0x1 03  ms1 
s,  =  29.9  psu 
H  =  1.0  m 
h  =  30.0  m 


u 


10 


=  5.0  ms1 
v10  =  5.0  ms"1 
Ta  =  -  30.0°  C 
T2  =  -  1.4°  C 
S2  =  30.8  psu 
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Run  2: 


INITIAL  CONDITIONS 


u,    =0.1  ms"1 
s,  =9.0  psu 
v,    =0.1  ms'1 
u^  =  3.0x10 3  ms'1 
vw  =  3.0x1 03  ms'1 
s1  =  29.9  psu 
H  =  1.0  m 
h  ■  30.0  m 


u10  =  10.0  ms1 
v10  =  10.0  ms 


T.  =  -  30.0°  C 
T2  =  -1.4°C 
S2  =  30.8  psu 
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Run  3: 


INITIAL  CONDITIONS  : 


u,    =0.1  ms"1 
s,  =  9.0  psu 
v,    =0.1  ms'1 
l^  =  3.0x10  3  ms"1 
vw  =  3.0x1 03  ms"1 
s1  =  29.9  psu 
H  =  1.0  m 
h  =  30.0  m 


u10  =  5.0  ms'1 
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Run  4: 


INITIAL  CONDITIONS 


u,    =0.1  ms 

s,  =  9.0  psu 
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Run  5: 


INITIAL  CONDITIONS 
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Run  6: 


INITIAL  CONDITIONS 


u,    =0.1  ms"1 
s,   ■  9.0  psu 
v,    =0.1  ms'1 
i^  =  3.0x10  3  ms'1 
vw  =  3.0x1 03  ms1 
s,  =  29.9  psu 
H  =  1.0  m 
h  =  15.0  m 


u10  =  5.0  ms"1 
v10  =  5.0  ms'1 
TB  =  -  30.0°  C 
T2  =  -1.4°C 
^  =  30.8  psu 
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Run  7: 


INITIAL  CONDITIONS 


u, 


=  0.1  ms1 


*i 

=  9.0  psu 

v, 

=  0.1  ms1 

"w 

=  3.0x10^  ms1 

vw 

=  3.0x10^  ms1 

*1 

=  30.9  psu 

H 

=  1.0  m 

h' 

=  45.0  m 

u 
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=  5.0  ms1 
v10  =  5.0  ms"1 
T.  =  -  30.0°  C 
T2  =  -1.4°C 
^  =  31.5  psu 
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Run  8: 


u, 


=  0.1  ms1 


*l 

=  9.0  psu 

v, 

=  0.1  ms1 

"w 

=  3.0x103  ms1 

Vw 

=  3.0x10  3  ms1 

Si 

=  29.9  psu 

H 

=  1.0  m 

h  = 

=  30.0  m 

INITIAL  CONDITIONS 


"io  = 

5.0  ms 

Vio   = 

5.0  ms 1 

T.  = 

-  30.0°  C 

T2  = 

-1.4°C 

Sj  =  31 .25  psu 


IICWIK  IMCUtt  IUI«1HM 


It* 


III 


til 


151 


HI 


„_£. 


It — 


± 


?l 


II 


41 


SI 


tan 

OVTCOIK  10IHAVE   lAtUTIM 


MM 


71 


COIPVCTIVE   *€»!   IIVX 


Nil 


fERtlltl  «»T  fl«X 


•M 


1.11 
t.M 
l.tl 
1.11 
I.II 
I. It 
1.71 


__4_ 

■     : 

_| — 

_H 

I\ 

j — 

f-.-^-S 

4. 

— j 

f_ 

: 

:     - 

-     \  ^S* 

.[_     __... 

-    -:■■■        ..|. ............  ■■ 

_  ..  _j_.    _    ._. 

1                           .      "     . 

' "  "1  "  - - ■ 

.    "•' 

?• 


w» 


BCfM   HU1   HUI 


3            -                         — ! ' ' ■ ' ! ' ' ' ' ! ' ' ' ~"T '        *        ! 

\ ^1  [ 

_-— — 4 ^^                                         = 

— 4— -~ — " 1               \               \               \  '- 

: 

'•                                             !                                              !                                              I                                              >       " 

; 

:                              :                              :                              1                              1    - 

L 

i            i            1            1            1  = 

11 


u 


it 


Dim 


72 


cimiMKii  mocin 


DtJI 


•I.  II    — 


1.14 


TfMffHTWl   MIFHEtCf    (TIT?) 

-r 


-1.1) 


I. II    — 


■l.fl  I — h 


•*!• 


SUIIITT   filMEHilCf    (fl(l) 


hm 


73 


UMMCf   TEtVEUTME 


ICE   THICTREIS 


I 

l 

1 

!            t            !  ^-- — T 

■ 

'-- 

: 

■ 

': 

™ 

| 

■ 

,  .1 , 

MIXED  UTER  MHH 


41 


74 


Run  9: 


INITIAL  CONDITIONS 


u,    =0.1  ms'1 
s,  =  9.0  psu 
v,    =0.1  ms'1 
u^  =  3.0x10^  ms1 
vw  =  3.0x10  3  ms'1 
s1  =  29.9  psu 
H  =  1.0  m 
h  ■  30.0  m 


"io 

=  5.0  ms'1 

Vm 

=  5.0  ms 1 

T. 

=  -  30.0°  C 

T? 

=  -1.4°C 

S2  = 

=  30.35  psu 
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Run  10: 


INITIAL  CONDITIONS 


u,    =0.1  ms 
s,  =  9.0  psu 
v,    =0.1  ms'1 
=  3.0x103  ms 1 


=  3.0x10  3  ms1 


s1  =  29.9  psu 
H  =  1.0m 
h  =  30.0  m 


u10  =  5.0  ms"1 
v10  =  5.0  ms"1 
T.  =  -  30.0°  C 
T2  =  -1.2°C 
Sj  =  30.8  psu 
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Run  11: 


INITIAL  CONDITIONS 


u. 

=  0.1  ms1 

s, 

=  9.0  psu 

v, 

=  0.1  ms'1 

"w 

=  3.0x10  3  ms1 

K 

=  3.0x1 03ms1 

&1 

=  30.35  psu 

H 

=  1.0  m 

h  = 

=  30.0  m 

u10  =  5.0  ms'1 
v10  =  5.0  ms"1 


T.  =  -  30.0°  C 
T2--1.4°C 
&2  -  30.8  psu 
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